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Abstract 

The widths for tttt decay of the L = 1 charm mesons are calculated by describing the pion coupling 
to light constituents quarks by the lowest order chiral interaction. The wavefunctions of the charm 
mesons are obtained as solutions to the covariant Blankenbecler-Sugar equation. These solutions 
correspond to an interaction Hamiltonian modeled as the sum of a linear scalar confining and a 
screened one-gluon exchange (OGE) interaction. This interaction induces a two-quark contribution 
to the amplitude for two-pion decay, which is found to interfere destructively with the single quark 
amplitude. For the currently known L — 1 D mesons, the total tttv decay widths are found to 
be ~ 1 MeV for the Di(2420) and 3 MeV for the 1)5(2460) if the axial coupling of the constituent 
quark is taken to be = 1. The as yet undiscovered spin singlet Di state is predicted to have a 
larger width of 7 - 10 MeV for tttt decay. 
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1 Introduction 



The pion decay widths of excited charm mesons, the D mesons, which are formed of one hght flavor quark 
(antiquark) and a charm antiquark (quark) provide the best observables for determining the couphng of 
pions to hght constituent quarks, as the pions do not couple to the charm quarks. The couphng strength 
of the pion to the hght flavor constituent quark is determined by the axial coupling constant of the 
quark. The theoretically indicated values for ||, |[ IJ], which fall in the range 0.75 to 1.0, have 

been shown to be consistent with extant empirical information on the single pion decay widths of the 
ground state and excited charm mesons 01. The uncertainty margin derives from the fact that only 
upper limits have been set experimentally M on the widths of the D* mesons, while the total decay 
widths of the orbitally excited L = 1 Z3i(2420) and 152(2460) charm meson states are known only within 
a very wide uncertainty range. As the latter lie well above the threshold, not only for single pion but 
also for two pion decay, it would be particularly instructive to obtain theoretical predictions for as well 
as empirical information on the branching ratios for the latter decay modes. 

We here report a calculation of the two-pion decay widths of the excited L — 1 charm meson states, 
by extending a similar earlier calculation of the widths of their single pion decay modes |^ . The charm 
mesons consisting of a light quark (or antiquark) q and a heavy antiquark (or quark) Q are here treated 
as relativistic two-particle systems described by wavefunctions which have been obtained as solutions to 
the covariant Blankenbecler-Sugar equation || in rcf. The interaction between the q and the Q (or 
the q and the Q) is described as a combined linear confining and screened one-gluon exchange interaction, 
with parameters determined so as to obtain a satisfactory prediction of the empirically known part of 
the D meson spectrum JTo[ . 

The model for emission of two pions from the light constituent quarks employed is the conventional 
chiral pion-quark pseudovector coupling model, which includes a Weinberg- Tomozawa type term for 
constituent quarks. The only parameter in this amplitude is the axial coupling of the light constituent 
quarks , as the quark mass is fixed by the Hamiltonian model used to determine the D meson spectrum 
and wavefunctions. In addition to the single quark amplitude for two-pion emission, we also consider 
the interaction current contribution to the vector current of the qQ system, which is associated with 
intermediate antiquark terms. This exchange current contribution to the Weinberg- Tomozawa interaction 
is found to interfere destructively with the single quark amplitude and to bring about a reduction of the 
calculated width for the two pion decay of the L = 1 D mesons by about 25%. 

The amplitude for two-pion decay is found to be mainly spin-independent, which is in accordance with 
the current empirical information on the tttt decays of the analogous strange K mesons. The tttt decay 
width of the meson is predicted to be about 3 MeV. For the Di meson, the corresponding number 
is found to be 1 MeV. Estimates are also given for the tttt decay widths of the hitherto undiscovered D 
meson states with L ~ 1, i.e. the spin triplet Dq and the spin singlet U*. The masses of these two states 
are predicted in ref. to be 2340 and 2390 MeV respectively. The tttt decay width of the Z3g is here 
found to be small (of the order one tenth of an MeV, whereas that of the Z?* is predicted to be about 
7 MeV. This variation is a natural consequence of the rather small phase space available for these decays, 
combined with the significant spin-spin and spin-orbit splittings between various D meson states. 

This paper is divided into 4 sections. In section 2 the operator for two-pion emission from a single 
quark is described along with a derivation of the corresponding decay widths. In section 3 the exchange 
current contribution to the two-pion decay width is calculated. Section 4 contains a concluding discussion. 
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2 The Decay Rate for tttt decay 



2.1 Single Quark Amplitude for Two-Pion Decay 

The emission of two pions from a D meson may be described by the pseudovector Lagrangian, which 
constitutes the lowest order chiral couphng for pions to constituent quarks: 

g1 _ 

^ = ^of"^? ^57m df, 4>7r ■ T^jjq. (1) 

Here g^^ denotes the axial coupling constant of pions to light constituent quarks, and is the pion decay 
constant, the empirical value of which is 93 MeV. The axial coupling constant is conventionally taken 
to be equal to, or somewhat less than, unity This coupling gives rise to Born term and crossed 

Born amplitudes of conventional form. Fig. |l| (a,b), for the emission of two pions from an interacting 
constituent quark. 

The chiral model for the two-pion decay amplitude is completed by the Weinberg- Tomozawa (WT) 
interaction, which is described by the Lagrangian 

1 — -* -* 

CwT = -JT^'ipqlp.f ■ (j)^ X df,(l)ni^q. (2) 

J TT 

This interaction leads to the contact amplitude that is described by the diagram c) in Fig. |. The Born 
amplitude (Fig. ^a) may be expressed as 

/ *^ \ ^ 1 

= -M ^77" ) i'qlbl-h : 'y57 ■ kalpqTbTa- (3) 

V 2/7r / 1-Pa- 

Here denotes the mass of the light constituent quark, which in rcf. [0 was obtained as 450 MeV, 
as this value led to an optimal description of the charm meson spectrum within the model considered. 
Similarly the crossed Born amplitude (Fig. |l|b) takes the form 

TcB = -i \ 777" ) Tpq 757 • ka : 757 ' h ^qTaU. (4) 

The general isospin decomposition of the two-pion emission amplitude for constituent quarks is, in 
analogy with that for nuclcons, 

T^SabT+ + ^[n,Ta]T~. (5) 

The general invariant expression for for the amplitudes and is in turn 

T± = u{p') (v4± - ij ■ QB^) u{p). (6) 

In these expressions, the four- vector Q denotes the combination Q — {kb — ka)/^, where ka and kb are 
the four-momenta of the emitted pions. The commutator of the SU(2) generators may be expressed as 
[Tb,Ta] — 2ieabcTc- In this notation the Born, crossed Born and Weinberg- Tomozawa amplitudes are, 
respectively 

= ' (||) + + 4"^'^2^) {^"a + I [n, r.]) , (7) 
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Figure 1: Feynman diagrams that correspond to the interaction Lagrangians of eqs. and (|^). The 
diagrams a) and b) describe the Born and crossed Born ampUtudes and the diagram c) describes the 
contact vertex that is associated with the Weinberg- Tomozawa Lagrangian, eq. (^. In these diagrams 

Pa ^ p ~ ka, Pb = P - h and ka + kb ^ p - p' . 

Comparison of these amphtudes with eq. yields the desired expressions for the sub-ampHtudes 
and , which are 
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Here the identities = — s and ~ ~u, where s and u are the invariant Mandelstam variables, have 
been used. These results for the A and B amplitudes are formally equivalent to the corresponding results 
for the two-pion emission amplitude for nucleons Note that in eq. (|l^), the contribution from the 
Weinberg- Tomozawa interaction tends to cancel the constant term in the B~ amplitude that arises from 
the Born terms. If the axial coupling constant g\ is taken to equal 1, then this cancellation is exact. 



2.2 The Width for Two-Pion Decay 

When excited D mesons decay by emission of pions, it is only the light constituent quark component 
that couples to the pions. Although pion decay from a single non-interacting quark is kincmatically 
impossible, the pion emission amplitude nevertheless involves the charm antiquark indirectly through 
the bound state wave function in the matrix element of the single quark pion emission operator. The 
role of two quark emission operators that involve the interaction between the light quark and charm 
antiquark, with pion emission from the former, will be considered below in section ^. 

The general expression for the two pion decay width of an excited D meson may be written in the 
form: 

r = l^'f I W + + - fl). (14) 

J (27r)'^ (27r)-* (27r)'^ ^uJaOJb 

Here ka and fcf, denote the four-momenta of the two emitted pions, Pi and Pf denote those of the initial 
and final state D mesons while uoa and lo^ denote the energies of the emitted pions. Since the constituents 
of the D mesons form bound states, their normalization factors are included in the spinors u{p') and 
u{p) in eq. (||). In the laboratory frame P° = Mi. By introducing the variables Q = {ki, — ka)/2 and 
q = kfy + ka, the decay width expression may be rewritten as 



si^sff^Mj+Ua + ujb-M.y (15) 



Here the energy factors are defined as Wa = y rn^ + (9/2 — QY and uj^ = \J ml + {q/2 + QY respectively. 
The remaining delta function may be used to fix the variable Q, so that finally the expression for the 
differential width becomes 



d^,- H2^Yj, '^'"^ l/"c.,(<z,z)(g, + f)+c.,(,,z)(Q,-f) ■ ^^^^ 

Here it is understood that in order to obtain the total widths for two-pion decay, eq. (^) is to be inte- 
grated over fig yielding an additional factor An. In eq. (p^, the variable z is defined hy Q-q = Qqz. With 
this notation the pion energies oja and ujb are given by the expressions uJa — + + (?^/4 — Qfqz 



and ujb = Jyri^ + + (7^/4 + Qfqz, where the fixed variable Q/ is given by 



, _ [Ef - M,)4 - (Ami + <f){Ef ~ h'h? 

4{Ef-M^Y-4q^z^ ' ^ ' 



Here Ef denotes the energy of the final state D meson and is given hy Ef — ^ q^ + M'j. In all these 
expressions, ttItt denotes the pion mass. Since different charge states are not considered in this paper, 
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an average pion mass of 137 MeV has been used throughout the calculations. The cutoff momentum 
Qf corresponds to the maximal momentum of any one of the final state particles, e.g. the final state D 
meson. Thus g/ corresponds to the q-value of a decay of the form D' DX, where D' and D are the 
appropriate D meson states and X is a particle with mass Mx = 2m7r. The appropriate values for qf 
are listed along with the employed D meson masses in Table § 

The decomposition in eq. (||) is convenient with respect to the summation over the isospins of the 
pions, leading to 

|T/,p = 3T+^T+ + 6T-^T-. (18) 

isospin 

For calculational purposes, it is useful to split eq. (||) into spin independent and spin-dependent parts as 

=a^ (3^. (19) 

Here, the amplitude a will contain, in addition to the spin independent amplitude A, the spin-independent 
contributions that arise from the 7 • Q term in eq. (|^) , while f3 will contain all the spin-dependent terms 
contained in the expression eq. (^. Consider first the non-relativistic limit. In that limit, the expressions 
for a and /3 become 



(20) 



q X 



2m„ 



(21) 



where Qo is defined as {ujb — '^a)/2 and P is defined as P = (p^ + Pq) /2. Because two pion decays of 
charm mesons with L = 1 to the ground states with L = are considered here, the largest contribution 
is expected to arise from the P dependent term in the spin independent amplitude a^. 

Evaluation of the expression (|8|) requires spin sums for both the spin independent and spin-dependent 
terms in eq. ( p^ ) . The spin independent part of the amplitude requires evaluation of the spin sums 



1 



2J+1 



J2 {LSJM\a* 



M=-J 



1 



^]a\LSJM), 



(22) 



1 



2J+1 



J2 {LSJM\a* 



M=-J 



a \LSJM) 



(23) 



for decay to the spin singlet D and spin triplet D* mesons respectively. For two-pion decay of the spin 
triplet Z?i(2420) and 132(2460) mesons the spin independent amplitude gives no contribution for decay 
to the ground state D meson (p^). The contribution given by the expressions above is \a\^ to the decay 
rate to the spin triplet D* meson (p^). 

The corresponding spin sums for the (iq-P dependent term in eq. (Q^ are somewhat more complicated. 
For initial states with L = 1, S" = 1 they may be expressed as in ref. [H, giving 
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y (iijA/i^- 



f3' 5i2(/3) 



for decay to spin singlet (D) final states, and 

^ .7 

2J 



6 



IIJM), 



1 



2/32 , 5i2(/3) 



UJM) , 



(24) 



(25) 



for decay to spin triplet (D*) final states. In eqs. ( p^ and (p5|), S'i2(/3) is defined as the tensor operator 
Si2{P) = 3tTg • /3(Tq ■ (3 ~ Uq ■ ctq. Givcn these spin sums, and the state vectors for the spin 1 P-state D 



mesons, 



\isjM) = ^ (ms\jM) ^^^Yu{f) \is) + {im\jM) ^^yi/(f) |oo) 



(26) 



where \\s) denotes a spin triplet state with Sz — -s, the following spin summed squared amplitudes for 
the two-pion decays of the spin triplet D2, Di and Dq mesons are obtained: 



\T\ 



3 

3 [I 
3 
0. 



a+f + 2\a-\'' + ^{\0+\'' + 2\f3 
'+\^ + 2\a-\^ + 0+\^ + 20-\^ 

+p+2Kn, 



1(1/3- 



2ir 



\f3+\' + 2\r\' 



(27) 

(28) 
(29) 

(30) 

(31) 
(32) 



From the above results, it follows that for the tttt decays of the spin triplet Di meson, there is no 
contribution for two-pion decay to the D states. For the spin singlet DI meson, similar expressions hold, 
with the modification that the spin-independent amplitudes now only contribute to decay to D states. 
For the spin-dependent terms, one has to consider the following additional spin sums 



i ^ (lOlAfI 



5i2(/3) 



llOlAf), 



Af=-1 

for decay to spin singlet D final states, and 

1 

3 



i J2 (lOlMl/? 
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llOlAf) , 



(33) 



(34) 



for decay to spin triplet D* final states. Application of these spin sums yields then the desired expressions 
for the spin summed squared amplitudes for two-pion decay of the spin singlet meson: 
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\T?r 



DT^D 



= 3 
- 3[|a 



l/3+r + 2|/3-|^ 



(35) 
(36) 



In the non-rclativistic limit, the foUowing two radial matrix elements are required for the numerical 
evaluation of the above amplitudes: 



Mo 
Ml 



3m. 



dr 



q Jo 



Uo(r)Mi(r) - uo{r)u[{r) 



Uo{r)ui{r) 



Jo 



(f) 



I dr uo{r)ui{r) ji (^^) ■ 



(37) 
(38) 



Here uq and ui denote the reduced radial wavefunctions that arc obtained by numerical solution of the 
Blankcnbecler-Sugar equation in ref. for states with L = 0, 1 respectively. 

In the following calculation, the non-local combinations of the kinematic variables, which appear in the 
denominators in the Born term amplitudes and B~ , have been approximated by their corresponding 
expectation values. For calculational convenience, the approximation Eg ~ rriq has likewise been made. 
The kinematical variables Q and q have, in addition to the energies of the emitted pions, been treated 
without approximation. The resulting expressions are then obtained as 



B- 



-4m 



, q \ 2 
2 / 9 A 



-4m 



, q \ 2 
2 / 9a 



2U 



P-q)-2(P 



Q ■ q- 2mqUJa 



P-q) + 2(P 



— + Q ■ q~ 2mqUJb 



P-q)-2{P-Q 



P-q)+2(P 



— + Q ■ q - 2mqUJb 



- Q ■ q- 2mqUJa 



T q 2 



(39) 



(40) 



The contribution to the decay rate from the spin independent a"*" amplitude in ( |l9|) may then be expressed 
as 



\T\l^ = 'i\a+\' = 'i 



A-Mi+B+i^^M,-^Mo 



(41) 



for all decays to D* final states. Note that the matrix elements here are (/-dependent. Similarly, by using 
the B^ amplitude, the contribution to jTp from the term becomes 



\T\'=6\a- 



6 (B- 



— T — Ml ^A^o 



(42) 



8 



for all decays to D* final states. Here uja and cjf, are defined as for eq. (p^). In addition, the i?+ and B 
amplitudes also contribute through the spin dependent amplitudes and f3~ (|l^)- In the case of the 
decay mode D2 — > I?*7r7r, these contributions may be expressed as 



18, 



3 9^,(1 



and 



B- 



(43) 



5 4 mi 



B^ 



(44) 



The corresponding expressions for the remaining decay modes of the spin triplet L = 1 D mesons to which 
the spin-dependent amplitudes contribute may be inferred from eqs. ( p7[]3^ ): For the decay — *■ Dtt-k, 
the factors 3/5 on the r.h.s. of eqs. ( |4^ ) and ( ^ ) are to be replaced by 2/5, and for the decays Dq Din: 
and Di — > D*Tnr, those factors equal 1. That is also the case for the decay of the spin singlet D* meson 
to D*7r7r (|36|). The total decay widths are obtained by adding the contributions from eqs. ( p| , ^ 
and (0), and inserting them into eq. (p^. In the numerical evaluation of the expressions ( ^9| ) and (4^, 
the terms P ■ q and P ■ Q in the denominators have been approximated by their respective expectation 
values: 



P-q 



2{P-Q 



QfzrriqMo. 



(45) 
(46) 



Note that despite the appearance of the matrix element Mq in the above expressions, no non-relativistic 
approximation is implied. 

As the velocities of the confined quarks and antiquarks in the D mesons are close to that of light, the 
non-relativistic approximation for the amplitudes is not reliable. Therefore, eqs. (|20| ) and ( ^l|) should be 
replaced by the analogous unapproximated forms 



2E' 



PQ 



2E 



A± ( 1- 
1 



{E' + mq){E + mq) 
Q-q( 1 1 



E + niq E' + mqJ 2 \E + ruq E' + niq 




P^-q^/A 



(E' + mq){E + mq) 



(47) 



and 



0± ^ jE' + mqjE + mq 



2E' 
qx Q 



2E 



{E + mq){E' + rUq) 
1 



qx P + B^ 



PxQ 



E + rtiq E' + niq 



1 



1 



E + WLq E' + IJlq 



(48) 
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Thus, relativistic counterparts to eqs. (^^^^ are required. As a consequence, eqs. and (H) should 
be replaced by the expressions 



and 



ITI 



6(B- 



Ml' 



Here the relativistic matrix elements A^™' arc defined as 



(49) 



(50) 



1 pOG nOO nOQ nl 

M'^'l = - dr'r'uoir') / drrui{r) / dP dvfBs{P,v) 
Jo Jo Jo J~l 



q/A + Pv 



^P^ + qyi6 + Pqv/2 V 2E' 



E' + TTig E + nig 



2E 



{E' + mq){E + mq) 



(51) 



where the last term in eq. ( p7|) has been neglected because of its smallness. Similarly, the relativistic 
matrix element A^ff' is obtained as 



-Mo 



rcl 



-f pOO pOC pQO pi 

— / dr'r' / drr dP P^ dvfBs{P,v) 

•JTT Jo Jo Jo J-1 



2E' 



E' + nig E + niq 



2E \E + niq E' + rriq 



u'air')ui{r)-uoir')u[{r)^2 



Uoir')ui{r) 



(52) 



In the non-relativistic limit, these matrix elements reduce to the forms of eqs. ( |38|) and (^7|) respectively. 
In the above equations, the factor /bs(-P, i^) arises in the reduction of the amplitude from the form 
appropriate to the Bethe-Salpeter equation to that of the Blankenbecler-Sugar equation, and is defined 
as 



fBs{P,v) = 



M, 



Q 



(53) 



^{E + Ec)iE' +E'^) 

Here Mq denotes the mass of the heavy (anti)quark. The energy factors E,E' and Ec,E'^ are defined as 



E = J ml + P2 + pqy + q2 j^^ e' = J m? + P^ - Pqv + q^ /4, 



(54) 



for the light quark g, and 
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(55) 



for the heavy (anti)quark Q. In addition, eqs. (43) and ( [t4[ ) should also be replaced by analogous 
relativistic forms. These may be obtained as 



\T\l^=&{B-f 



H 2 



for the contribution from /? , and 



■ -Mm 



/32 



for the corresponding one from Here the matrix elements are defined as 



(56) 



(57) 



dr'r'uo{r') / drrui{r) / dP / dvfBsiP,v) 



q/4: + Pv 



^P^ + q^/16 + Pqv/2 V 2E' 



E' + rriq E + rUg 



2E \E' + mq E + niq 



(58) 



and 



dr'r'uoir') / drrui(r) / dP P^ / Vl - i^Vbs(^', ^') 
^0 ^0 



q/A + P% 



E' + mq E + mg 



^/P'^ + q'^/16 + Pqv/2 V 2E' 



JO r'Jp^ + ^-^]n{r^|P^ 
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1 



2E {E' + mg){E + mg) 




(59) 



With exception of a factor , the matrix element (|5^) reduces to the form ( |38D in the non-relativistic 
limit. In these expressions, the P x Q term in eq. ( ^8| ) has been dropped because of its smallness. Fur- 
thermore, it turns out that eq. ( p9[ ) is numerically much smaller than eq. (^). Overall, the contributions 
from the amplitudes a'^ and a~ are dominant, whereas those from /3+ and f3~ represent only small 
corrections. This feature is in accordance with the current experimental status for the analogous decays 
of the excited strange mesons, for which the decay mode i^l Ktttt has not yet been seen, although 
the decays of the K2 meson are otherwise well studied. In eqs. (^6|) and (|5^), it is again understood that 
the factor 3/5 is to be replaced as for eqs. (|4^) and ( ^ ) when different decay modes are considered. The 
constituent quark and D meson masses used are listed in Table |l|. The numerical results obtained by 
using the relativistic expressions are given in Table H For comparison the corresponding non-relativistic 
results are presented in Table ||. The non-relativistic treatment leads to overpredictions by a factor ~ 3. 
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Decay 


Initial state mass [Ali) 


Final state mass (M/) 


If 




2460 


2007 


327 


Dtttt 


2460 


1867 


462 


Dl D*TTTT 


2420 


2007 


282 


Dl D*TTTT 


2389 


2007 


244 


Dl Dtttt 


2389 


1867 


395 


D* D*TTTT 


2341 


2007 


177 


Dl Dtttt 


2341 


1867 


347 



Table 1: Initial and final state D meson masses and resulting maximal kinematically allowed momenta 
qf in MeV used in the calculations. The masses listed are to be viewed as averages over the different 
charge states. In the calculations, the values mq = 450 MeV and Mq = 1580 MeV, obtained by fits to 
the spectra in ref. [0, have been used for the constituent quark masses. 



Decay 


\T\l. 


1 la 


\T 


2 

P- 




\T 






Total {g\ = 1) 


Total {g\ = 0.87) 


Dl -> D*TTTT 


0.896 


2.864 


5.06- 


10- 


3 


1.17- 


10 


-3 


3.77 MeV 


2.39 MeV 


Dl -> Dtttt 






6.45- 


10- 


'2 


7.51 • 


10 


~3 


0.07 MeV 


0.05 MeV 


Dl D*TTTT 


0.367 


1.291 


2.33- 


10- 


3 


7.85 • 


10 


-4 


1.66 MeV 


1.05 MeV 


Dl D*TTTT 






6.54- 


10- 


4 


3.30- 


10 


-4 


~ MeV 


~ MeV 


Dl Dtttt 


2.749 


7.915 














10.7 MeV 


6.80 MeV 


D^ D*TTTT 


0.020 


0.098 














0.12 MeV 


0.07 MeV 


Dl Dtttt 






1.43- 


10- 


2 


2.86- 


10 


-3 


0.02 MeV 


0.01 MeV 



Table 2: Numerical results using the single-quark amplitudes and the relativistic matrix elements for the 
two-pion decay widths of the spin triplet D^, Di, D^ and the spin singlet Dl mesons. The individual 
results from the spin independent (a^) and spin dependent amplitudes (/3^) are shown for 5^ = 1, and 
the total TTvr decay widths for both — 1 and — 0.87. 

The results for the two-pion decay width of the positive parity charm mesons with L = 1 shown in 
Table || reveal a strong sensitivity to the strength of the pion coupling to the light constituent quarks 
as measured by the axial coupling constant of the constituent quarks. This is of course as expected, as 
the calculated decay widths are proportional to the 4th power of g^. Also, because of the large (140 
MeV) splitting between the D and D* states and the significant spin-orbit splittings of the L = 1 charm 
mesons, it turns out that the values of the maximal kinematically allowed momentum transfer qf as 
given in Table || show marked variation. Consequently, some decays are kinematically favored, while 
others, in particular D^ — *■ D*tttt, are strongly inhibited by the small phase space available. Hence 
the two-pion decay widths of the L ~ 1 D mesons are very sensitive to the spin-orbit structure of 
the quark-antiquark interaction. The same conclusion was also reached in ref. ||^ concerning the single 
pion decays of the L ^ 1 D mesons. In the absence of empirical information and definite QCD lattice 
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calculations ^ |T^, |l^ , the energies of the hitherto undiscovered Dq and mesons have here been 
taken to equal those found in the calculation of ref . |^ . 

The predicted total decay widths of the L — 1 D mesons in the single quark approximation may 
be obtained by adding the calculated two-pion decay widths in Table ^ to those for single pion decay 
obtained in ref. With g\ = l the total calculated tttt decay width of the 1)2(2460) is 3.8 MeV and 
that of the Di(2420) is 1.7 MeV. If these values are added to the corresponding calculated values for 
single pion decay obtained in ref. |6|, the total calculated width of the _D2(2460) comes to 19.5 MeV, 
which is well within the uncertainty margin of the empirical value 25^^ MeV for the total decay width [Q . 
In the case of the i3i(2420) meson the total calculated width for single and two-pion decay comes to 
15.3 MeV, which is close to the empirical uncertainty margin of the total decay width IS.Qlg'g MeV 0. 
Reduction of the value for g\ to 0.87 brings the calculated values for the total width for it and tttt decay a 
bit below the empirical values for the total widths. Likewise, employment of the two-quark contribution 
considered in the next section has the effect of reducing the calculated tttt widths. The final results for 
the TTTT widths are listed in Table ||, along with predictions for the total widths of the L = 1 D mesons. 



Decay 


\T\l. 


\T?- 

1 la 


\T\l- 


\T 


2 

13+ 




Total {g\ = 1) 


Total {g\ = 0.87) 


D* ^ D*7:tt 


2.638 


9.934 


1.36- 10~ 


-2 


3.11 - 


10 


-4 


12.6 MeV 


8.01 MeV 


Dl Dtttt 






0.178 




5.77- 


10 


~3 


0.18 MeV 


0.12 MeV 


Dl D*TTTT 


1.083 


4.539 


6.21- 10" 


-3 


1.16- 


10 


~4 


5.63 MeV 


3.56 MeV 


Dl D*TTTT 






1.74- 10" 


-3 


2.60- 


10 


-5 


~ MeV 


~ MeV 


Dl Dtttt 


8.048 


26.97 












35.0 MeV 


22.4 MeV 


D* D*TTTT 


0.059 


0.357 












0.42 MeV 


0.26 MeV 


Dl Dtttt 






3.84-10- 


2 


9.51 - 


10 


-4 


0.04 MeV 


0.03 MeV 



Table 3: Numerical results for the two-pion decay widths of the spin triplet D2^ Di, Dq and the spin 
singlet Dl mesons using the non-relativistic approximation to the single-quark amplitudes. The individ- 
ual results from and (3^ are shown for g\ = \, and the total decay widths for both g\ = \ and g\ 
= 0.87. 



This apparent slight underprediction of the total decay widths may actually be a desirable situation 
in view of the fact that it is kinematically possible for other decay modes, mainly tttttt and ry-meson decay, 
to contribute to the total widths of these mesons, even though there is very little phase space available 
for such decays. However, in view of the large systematical uncertainties involved in the experimental 
determination of the total widths of the Z?2 and Di mesons, the uncertainty margins quoted by ref. 
may be on the narrow side. 
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3 Two-Quark Contribution to the Decay Width 



It is instructive, for the purpose of determining the two-quark contribution to the decay width, to rewrite 
the Weinberg- Tomozawa Lagrangian, eq. (^) in the form of a current-current couphng: 

>CwT = -TTJ^M • ^ 9^$^. (60) 

Here = itpqjp^ipqT is the isovector current of the Hght constituent quark and (/)^ x S^^^r is the 
current of the two-pion system. Given this expression it becomes very natural to describe the irreducible 
two-quark contribution to the two-pion production operator by means of two-quark interaction current 
contributions to the isovector current Vfj,. 




Figure 2: Two-quark contributions to the pion production amplitude associated with the Weinberg- 
Tomozawa Lagrangian. The diagrams shown correspond to both time orderings of the two-quark contri- 
bution from the scalar confining interaction. 



The most significant interaction current contribution to the isovector current of the charm mesons is 
the "pair" current that is associated with the scalar confining interaction between the constituent quarks. 
This current, which is illustrated by the Feynman diagrams in Fig. |^, represents a relativistic correction 
term, which arises in the elimination of the small components of the quark spinors. In the non-relativistic 
approximation it can be described as a renormalization of the isovector current of the light constituent 
quark by the confining interaction, either in terms of interaction correction to the constituent quark mass 
or as an explicit additional term Q . 

In the non-relativistic limit, the isovector current — (iVo,V) of the light constituent quark takes 
the form 

V = -^{p'q+Pq-iaqXq)T, (61) 
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where is the hght constituent quark mass and Pq and p ^ are the initial and final quark momenta 
respectively. The corresponding expression for the "pair" current that is associated with the confining 
interaction is then obtained as 

V, = -Y^V. (62) 
rriq 

Here Vc{k) is the (formal) Fourier transform of the confining interaction, and the momentum transfer k 
is defined as the difference between the final and initial momenta of the charm (anti)quark: k — p'q —pq. 
Within the non-relativistic approximation it now becomes a straightforward task to take this two-quark 
contribution to the two-pion production operator [|l^ into account. The most significant contribution 
may be included by modifying that part of the amplitude in the expression for a" (eq. (|2^)), which 
arises from the Weinberg- Tomozawa Lagrangian, according to the following prescription: 



(63) 



Here Vc{r) is the scalar confining interaction, which in the wavefunction model of ref. |^ has the form 
Vc{r) ~ cr — b, with c = 1120 MeV/fm and b = 320 MeV. This replacement should only be made for 
terms that are associated with the spatial current coupling, i.e. for those that contain the momentum 
transfer Q. A similar modification should also be made in the expression for (3~ (eq. (|2l|)), but as that 
contribution is already very small compared to those from the spin-independent a terms, see Table ||, it 
will not be considered in this paper. 

The hyperfine interaction that is modeled as screened one-gluon exchange (OGE) between the quark 
and the antiquark also contains a central component which gives a contribution to the isovector exchange 
current. This contribution may be obtained as ||l0|] : 



Vgik) 



vn -* — ' 

-T7-{pq' +PQ +i^Q -X k) +iaq X k 



(64) 



In the above expression for the OGE exchange current, the Mg^ dependent term arises from the spatial 
current coupling 7^ • 7q term, and the second term from the charge coupling 7^7^ term in the the OGE 
interaction. Vg{k) denotes the form of the OGE interaction in momentum space, which is conveniently 
expressed as 

= (65) 

where a color factor of 4/3 has been included. In the interaction model used in ref. |p^ , the running 
coupling constant was taken to have the following screened form: 

"^(^') = lrin[(F + L^)/A^]' 

where the parameters mg and Ao were obtained as nig — 240 MeV and Aq = 280 MeV respectively. In 
case of the tttt decay of the D meson states with L = 1, the term with the sum of charm quark momenta 
in eq. (|6^), gives the largest contribution. As this term is inversely proportional to the mass of the 
charm quark, it is evident that the gluon exchange contribution will be smaller than the corresponding 
contribution from the exchange current (n2) that is associated with the confining interaction. This main 
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term of the isovector exchange current that is associated with the OGE interaction may be taken into 
account in the non-relativistic approximation by modifying eq. (p3h to read 



1 



1 



1 



ruq Mq 



(67) 



where Vg{r) is the central component of the one-gluon exchange potential. For a bare OGE interaction 
one would have Vg{r) = — 4as/3r. This simple expression is however significantly modified by the 
running coupling as given by cq. (|66|). This effect may be included as in ref. ||l0| , by expressing Vgir) as 



4 2 
3^ 



dkjoikr) Q!s(fc^)- 



(68) 



This form reduces to the static OGE potential given above, if the running coupling as is taken to be 
constant. Note that the relativistic corrections arising from the quark and antiquark spinors have not 
been included in the above expression, since it is more natural to place them in the expression for the 
matrix element instead. Moreover, comparison of these expressions with those obtained for the confining 
interaction shows that the net effect of the exchange current that is associated with OGE will have the 
same sign as that, which is associated with the confining interaction. 

In the expressions ( |62[ ) and (|4|) , one factor of 1 /niq represents the static limit of the propagator of 
the intermediate negative energy light constituent quark. In view of the results obtained in the previous 
section, this static limit is expected to give rise to a considerable overestimate of the two-quark amplitude 
contribution. A more realistic treatment may be obtained if the static propagator l/ruq is replaced by 
the symmetrized form 4/(2mg + E + E') as in ref. Q. Thus in the relativistic case, the appropriate 
replacement of the Weinberg- Tomozawa term would be the following extension of the replacement (pTj): 



1 



1 



1 - 



\E, + Mq)[E'^ + Mq) 



-Va 



\f+r'\ 



AE,E', 
4m„ 



\r + r'\ 



2m„ + E + E' 



1 - 



?V4 



1 



2mq+E + E' \E'^ + Mq E^ + Mq 



{E'^ + Mq){E, + Mq), 
(69) 



where the factors containing E'^ and E^ arise from the spinors of the heavy antiquark line in Fig. |[ To 
take the relativistic effects into account demands employment of the relativistic version of the matrix 
element A^q°' also in the case of the exchange current contribution. First of all, the contribution to 
the decay rate from the a~ term, eq. (^0|), should be modified to take into account the two-quark 
contributions to the Weinberg- Tomozawa interaction. This can be accomplished by replacing eq. (50) 
with 



\T\i- = 6 



B- 



V 2 



J_QfZ 



(MConf _ _^OGE) 



1 2 



(70) 



The appropriate forms for the relativistic two-quark matrix elements in the above equation can be 
obtained by multiplying the expression for A^q°' with the relativistic factors included in eq. ( |69|) above. In 
addition, the argument |r-|-r '|/2 of the potentials will be approximated by the expression 
The resulting forms for the two-quark matrix elements in eq. ([70|) are then obtained as 
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Conf 



1 
3^ 



dr 'I 



u^,(r')ui(r)-uo(r')u'iW-2 



Uo{r')ui{r) 



dPP' / ^dvfBs{P,v) 



2m, + E + E' 



{E + mg){E' + m,) 



I {E, + Mq){E'^ + Mq) / 



AEE' 



E + rtiq E' + 



+ MQ)iE, + MQ) 



Jo(r^/P- + -- — 



(71) 



for the confining interaction, and 



2mq + E + E' 




u'o{r')ui{r) - Uo{r')u[{r) - 2 
Arua 



uo{r')ui{r) 



1 



E, + Mq E'^ + Mq 



OO /• 1 

2 



(^ + mg)(£'^ + mq) 

' Jo\r 



+ 



E + niq E' + m. 




(72) 



for the OGE interaction. The numerical results that follow when eq. (70) is employed are displayed for 
the most important decay modes in Table 0. 







\T?- 




Total 


Decay 


Rel 


+Conf 


+OGE 


9\ = l 


g\ = 0.87 


^ D*mr 


2.864 


2.377 


2.144 


3.05 MeV 


1.82 MeV 


Di D*TnT 


1.291 


1.076 


0.974 


1.34 MeV 


0.80 MeV 


Dl Dtttt 


7.915 


6.535 


5.872 


8.62 MeV 


5.17 MeV 



Table 4: Numerical results for the most important tttt decay modes obtained upon employment of the two- 
quark contributions to the Weinberg- Tomozawa interaction. The modifications to the |Tp_ amplitude 
with g\ — \ are shown as follows: The column "Rel" gives the one-quark result for each decay rate 
(cf. Table |), in the column "+Conf" the contribution from has been added, and in "+OGE", 

the results that follow when both A^q°"^ and AIq"^^ are employed are given. The resulting total decay 
widths for tttt decay are also shown, for both g\ — I and f/^ = 0.87. 
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4 Discussion 



In the absence at the present time of experimental data on the two-pion decay widths of the D mesons, 
it is instructive to compare the calculated results with the empirical knowledge of the analogous pionic 
decay modes of the positive parity strange mesons. The if|(1430) state, which is the strange analog 
of the 1)2(2460) has a width of ~ 13 MeV for nn decay 0]. Given the smaller phase space and much 
narrower radial wavefunctions, the present calculated tttt decay width value of 3.1 MeV for the D2(2460) 
hence appears to be reasonable. Furthermore, in ref. 0, the branching ratio K2 K*tt is given as 
(24.7±1.5)% while that for K2 K*TnT is reported to be (13.4±2.2)%. Thus, the available experimental 
data indicates that for this particular decay mode, the width for tttt decay ought to be ~ 55% of the 
width for tt decay. If the current calculation is compared to that in ref. the ratio of D*-ktt 
to D2 D*TT is obtained as ^ 60% for g\~\. One may thus conclude that a width for D\ — > D*-k-k 
of 3 MeV, as obtained in the current calculation, is what would be expected by comparison with the 
corresponding decays of the strange K mesons. The present results also appear realistic in view of the 
fact that the empirical nonobservation of the decay mode — > Kixix indicates that the amplitude for 
two-pion decay is mainly of spin-independent character. 

The strange analog of the Z?i(2420) state is most likely a mixture of the if 1 (1270) and the A'i(1400) 
states. The two-pion decay modes of these two states are dominated by p meson decay. As meson 
decay is kinematically impossible for the 1)1(2420) as well as the other L = 1 D mesons, these decay 
modes have little bearing on the decays of the latter. Addition of the calculated decay widths for single- 
and two-pion decay gives predictions for the total widths of the L = 1 D meson states. These results are 
shown in Table H. 



D meson state 


TT width 


TTTT width 


Total g\ = 1 


Total g\ = 0.87 


Experiment 


D*2 

Dl 

D*o 
Dl 


15.7 MeV 

13.6 MeV 

27.7 MeV 
13.2 MeV 


3.05 MeV 
1.34 MeV 
- 0.1 MeV 
8.62 MeV 


18.8 MeV 

14.9 MeV 
28.8 MeV 
21.8 MeV 


13.7 MeV 

11.1 MeV 
21.0 MeV 

15.2 MeV 


25t^ MeV 
18.913^ MeV 



Table 5: Total decay widths for strong decay of the L — 1 D mesons that follow when the results of 
the present calculation of the tttt decay widths are added to those for single pion decay obtained in 
ref. The individual results for tt and tttt decay are shown for = 1. The experimental results, when 
available, have been taken from ref. 0. Note that the Di is a spin triplet state with J = 1 and 73* is 
the corresponding spin singlet state. 

As indicated by Table ^, the present calculation of the two pion decay widths of the L = 1 charm 
mesons completes the calculation of the corresponding single pion decay widths of ref. The calcu- 
lations show that, if the charm mesons are described as two-particle systems formed of one light flavor 
constituent quark and one charm antiquark, where only the light quark couples to pions, a fair descrip- 
tion of the presently known decay widths of these mesons may be obtained. Given the chiral form for 
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the coupling between the pions and the Hght constituent quarks, a better description of the total pionic 
decay widths obtains with the value g\ = 1 than with smaller values, even though the empirical data 
are still very crude at this time. If other decay modes, such as tttttt and rj decay should turn out to give 
appreciable contributions to the total decay widths of the positive parity charm mesons, a somewhat 
smaller value for may be favored. However, as tttttt decay has so far not been detected for the strange 
K2 meson Q and since the 77 decay width of that same state is empirically found to be very small, that 
possibility is apparently not supported by experiment at this time. From the calculations of ref. |Q, 
one may conclude that the width of the decay mode D2 Drj is probably less than 0.25 MeV. 

Furthermore there is the possibility of two-pion decay of the D2 and Di mesons through an inter- 



mediate D1 which is close to its mass shell, an effect which has been investigated by ref. |15 . There, 
this mechanism was found to contribute significantly to the two-pion decay widths, although the effect 
is very sensitive to the widths and spin-orbit splittings of the L = 1 D mesons and is thus very difficult 
to estimate. 

Finally it should be noted that the results are quite sensitive to the exact form of the amplitude, 
as given by eq. (^0|). Therefore, crude approximations have to be avoided in the case of this amplitude. 
If the B~ amplitude is treated non-relativistically, which implies dropping the nonlocal and Q dependent 
terms altogether, large overestimates of that contribution to the tttt decay widths will result. 
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